Shape Function Effects in B 
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Owing to the fact that ~ rnt Aqcd, the endpoint region of the charged lepton energy spectrum 
in the inclusive decay B Xctvi is affected by the Fermi motion of the initial-state b quark bound 
in the B meson. This effect is described in QCD by shape functions. Including the mass of the 
final-state quark, we find that a different set of operators as employed in Ref. ||5| is needed for a 
consistent matching, when incorporating the subleading contributions in B ^ Xqlvi for both q = u 
and q = c. In addition, we modify the usual twist expansion in such a way that it yields a description 
of the lepton energy spectrum which is not just valid in the endpoint region, but over the entire 
phase space. 



I. INTRODUCTION 

The theoretical machinery for the determination of 
|ycfc| from semileptonic _B-meson decays has reached a 
mature state over the last years. With the very precise 
data on the exclusive B D^*^£p£ channels as well as 
on the inclusive decays B — )■ Xcii'i from the B-meson 
factories the theoretical description has been improved 
so much that currently relative theoretical uncertainties 
for I Kb I of less than 2% are quoted \^^, yielding a total 
uncertainty of about 2% Q]. 

The nonperturbative corrections to the inclusive de- 
cays are parametrically of order l/m^ and hence are ex- 
pected to be smaller than in the exclusive channels, where 
the corrections are of order l/rri^. However, the inclusive 
rate depends on the 5-quark mass rub, which needs to be 
determined in a suitable scheme as precisely as possible. 

In addition, also the parameters of the heavy quark 
expansion are needed, such as Ai and A2 at order l/m^ 
and the corresponding parameters appearing at higher 
orders. These parameters are obtained experimentally 
by taking moments of the various inclusive distributions 
(such as lepton energy spectra or hadronic invariant mass 
distributions) , but the higher moments become more and 
more sensitive to higher-order corrections in 1/mf,, since 
the leading contribution to the nth moment is roughly of 
the order 

When determining the heavy quark parameters from 
the lepton energy spectrum, the higher moments become 
sensitive to the endpoint region of the spectrum. Us- 
ing the 1/mi, expansion the leading term is the partonic 
rate and still a smooth function, but already the first 
non-vanishing nonperturbative contribution exhibits an 
irregular behavior which is unphysical. The situation is 
in fact very similar to the one in 6 — > it transitions, where 
it is known that these singular contributions can be re- 
summed into a shape function. For heavy to light de- 



cays, such as the decays b uivg and b — > 57, the twist 
expansion, that is, the resummation of nonperturbative 
contributions, has been performed to the subleading level 
in the 1 /rrib expansion 0, 0, IE ■ 

It has already been noticed some time ago that the 
light-cone distribution of the B meson also has a signif- 
icant effect on the endpoint region of the lepton energy 
spectrum in i? ^ Xcii^i- This is due to the fact that 
numerically Tib Aqcd- The charm quark mass thus 

has to be counted as ^TObAqcD when performing the 
power counting. The endpoint region is known to be de- 
termined by the light-cone distribution of the initial state 
and has the width ^/rn}J^QCD, which happens to be of 
the same order as rric. Therefore, it is useful to consider 
the effects of the light-cone distribution of the B meson 
also in i? — » X^lvi. 

In the present paper we perform this analysis and com- 
pare with the standard expansion. As a by-product, we 
suggest a modified twist expansion which can be applied 
over the full phase space, incorporating the twist expan- 
sion in the endpoint region as well as the usual local 
expansion in the rest of phase space. 

Furthermore, performing the limit ttIc — > we discover 
an inconsistency in comparison with previous work j^. 
It turns out that at subleading order additional opera- 
tors are needed which are formally of leading order, but 
have coefficients of subleading order. Expanding into the 
usual local expansion we obtain the correct result for the 
terms of order l/m^, indicating that our result is consis- 
tent. As a consequence, compared to Ref. additional 
nonperturbative input in the form of new shape functions 
appears. 



A. Inclusive Decay Rate 

Neglecting the masses of the leptons the energy spec- 
trum of the charged lepton in the B rest frame is given 
via the optical theorem by 
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where y = 2Et/mi, denotes the rescaled lepton energy, 

_ Gl\V,,\''ml 
° 192^3 ' 

and q stands for either u or c. The "B expectation value" 
is defined as {0)b = {B\0\B)/2mB, where the QCD 
states are normahzed to 2ms- The operator Tg has the 
form 

= S^"^^^ / d'^e-'«-T[W/J(a:)M^,(0)], (2) 

where si and p£ are the lepton spin and momentum. The 
effective weak current is 

with Pl = (1 — 75)/2 and ue denoting the lepton spinor. 
Plugging this into Eq. Q yields 

T^^^(~-lmbL^BbUf^ with (3a) 

Laf3 = la^ilpPL, (3b) 

= -iJ d^xe'^P^-P''>--T[{Da"PLq)ix){qi^PLi^im]. 

(3c) 

Here, La/3 is the leptonic tensor and H^^ represents 
the inclusive g-quark-neutrino loop with the momentum 
transfer pb — pi. The 6-quark momentum pt contains a 
large part mi,v, where v is the velocity of the B meson. 
As usual, assuming that Q = mi,v — pi sets a perturba- 
tive scale, large compared to fc = — rribV ^ ©(Aqcd), 
we may perform an OPE of the transition operator Tq in 
powers of A = Aqcd /rrib U>MM 

B. Light-Cone Vectors and Power Counting 

Using appropriate powers of mb we may work with 
dimensionless variables, which are denoted by a hat, for 
example Q = Q/rub = v ~ pi. We also define A = 1 — y 
and p = m'l/ml. 

As already noted, the c-quark mass satisfies ^ 
mb Aqcd and thus p ^ C(A). The kinematic endpoint in 
the local OPE is given by A = p, the partonic endpoint. 
The endpoint region of the lepton energy spectrum is de- 
fined by A ^ C'(A): and it is well known [3 that in this 
region the local OPE breaks down. However, it has been 
shown that one niay still perform a light-cone or twist 
expansion El [Hill. 

To set up the light-cone expansion we use the velocity 
V and the lepton momentum pg to define a basis of two 
light-cone vectors, 

V = 7;i'n- + f^)i Pe^Een^^yn, (4a) 



satisfying = = and n ■ n = 2. The metric is 
decomposed accordingly 

rjt"" = ^n^'h" + ^n^'n" + yy'f^, (4b) 

and a generic four-momentum p can be written as 

= \p-n^ + \p+n'' + P^, (4c) 

where we defined 

P+ = n-p, p^ ^n-p, = Tj'^^Piy. (4d) 

The power counting 

0+-A^O(A), p^O{A), kr^O{K), 

(5) 

yields the standard twist expansion by expanding in pow- 
ers of A, taking also into account that p ^ C'(A). 

Note that this power counting becomes wrong for small 
lepton energies, since A becomes of 0(1). In this case the 
spectrum is described by the usual local OPE. However, 
as we shall discuss below, with a slight modification of 
the twist expansion it is possible to describe the spectrum 
also for small lepton energies. 

The relevant kinematic variable in the OPE is p — 
Pb — Pi — Q + k, the light-cone components of which are 

p+ = A + A;+, p_ = l + A;_, p± = k±. (6) 

Obviously all k^ dependence occurs in the combination 
A_|_ = A -I- A;4.. In the local OPE the complete k depen- 
dence is expanded. In particular, this produces terms of 
the form k^/A, which become large, of 0(1), near the 
endpoint. The twist expansion avoids these terms, since 
only the fc_ and k± dependences are expanded. 

Alternatively to Eq. (jsj , we may as well treat the com- 
plete dependence on A-|_ (and p) exactly. That is, we use 
the power counting 

fc_ ^ 0(A), fci ~ 0(A), (7) 

and only expand in powers of fc_ and k±_ from the very 
beginning. The order in A of a twist term given by Eq. (|7|) 
corresponds to the order of the first local term it contains, 
once the usual local OPE is performed. Eq. 10 defines 
what we call the modified twist expansion and holds for 
large and small lepton energies. Our results will therefore 
be valid over the whole lepton energy spectrum, provid- 
ing a correct interpolation between the usual twist and 
local expansions. 

Our modified expansion is a direct extension of the 
usual twist expansion. The latter does not expand A+ = 
A + k-^ because both A and fc+ are considered C'(A). 
Once part of the /c+ dependence is left unexpanded, it 
is consistent to keep the complete k+ dependence unex- 
panded, as it just means to keep correct small terms, 
which can now be treated exactly. This is similar to the 
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rriiiV+k — mhV+k\ /mjjV+k — ^ mi,v + k+L 




FIG. 1: Feynman diagrams for the zero- and one-gluon matrix elements (to leading order in ois 



local expansion, where p = m^/ml, although being nu- 
merically small, is always treated exactly, because it can 
be treated exactly. Once we exclude fc+ from the power 
counting and treat it exactly, there is no need to count 
A as 0{A) anymore. Thus we can treat it exactly much 
like p. This in turn extends the validity of our expansion 
down to low lepton energies (where numerically A is of 
order one which caused the breakdown of the usual twist 
expansion) . 

In other words, starting from the full expression (in- 
cluding all fc-components) both expansions expand the 
k_ and k±^ components. This is where we stop, while 
the usual expansion additionally neglects terms of sec- 
ond and higher twist order caused by A+, e.g., terms 
like A^. We keep all those "kinematic" twist terms, be- 
cause they are of 0(1) for low lepton energies. Note that 
we do not claim to include all second-order twist contri- 
butions, i.e., our spectrum is only correct to 0{A) in the 
endpoint region, but it is correct to 0{A^) for low lepton 
energies. 

Therefore, the difference to the usual twist expansion, 
defined by Eq. 0, is that the modified expansion auto- 
matically keeps all twist contributions that in the usual 
power counting are of higher order only because of addi- 
tional factors of A+ or p. These contributions are purely 
kinematic and do not require additional operators in the 
twist expansion, but appear only as higher-order terms 
in the OPE cocfhcients. Taking them into account yields 
a consistent result valid over the full region of lepton en- 
ergies. 



II. OPERATOR PRODUCT EXPANSION 

To keep things simple and to exhibit the structure of 
the OPE we will perform it in terms of QCD light-cone 
operators. Schematically, it has the form 



r, = y dc^^c(")(c^)o„(c^). 



(8) 



We expand to ©(A^) as defined by Eq. (j?)). With re- 
spect to the usual twist power counting this includes all 
contributions of leading and subleading twist, as well as 
some of second and even third order. In other words, we 
obtain the full coefficients of all operators with up to two 
covariant derivatives, which in particular retains all local 
terms up to local O(A^). To do so, we need to evaluate 
the zero- and one-gluon matrix elements of T^, depicted 
in Fig.^ which we do to leading order in a^. 



A. Leading Twist 

We first consider the case q = c and later take the limit 
nic — *■ 0. The zero-gluon matrix element of Tc yields the 
well-known result 



{b\Tc\b) = e{f-p)[i-0 



Ub7"PLUb. (9) 



In a second step the B expectation values of the C'„(w) 
are parametrized in terms of shape functions. 



Indices in round brackets are completely symmetrized, 

Extracting the terms of 0(1), 

{h\Tc\b) = 9{A+ - p)(c^°\-k+) + OiA)^UbrPLUb, 

we obtain the leading term in the OPE of Tc, 

Tc = J du;9{A^ - p) (ci") (^)Oo" {u) + 0(A)) , (10) 

where Aj^ = A — lo. The leading operator has the form 
0^{u)^b5CiD++Lo)^^PLb, (11) 
and its coefficient is 

Ci°)(c^) = (1 - 3i?2 + 2Rl)nc + (2 - 3i?„ + RDA^Uc, 

(12) 
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with iD = iD — rrihV and = p/Ai^. Note that in 
the modified twist expansion we keep the contributions 
proportional to fia , which would usually be considered as 
subleading twist due to the additional factor of A^^ . 

The B expectation value of Oq{uj) is given to leading 
order by the leading shape function, 

{0^{uj))b = ^v"{B^\bJ{iD+ + Lu)b,\B^) + 0{A) 

= V/H + 0(A), (13) 

where |i?oo) denotes the B meson state in the infinite- 
mass limit and the by are the static heavy quark fields. 

Together with Eqs. 1^ and H10() we find the lepton 
energy spectrum at leading order 

^ = 2Toy^9{y) J dc^0(A„ - p)Tp{A^)f{uj), (14a) 

where {R = p/A) 

rp(A) = (1 - 3i?2 + 2i?3) + (2 - 3i? + i?3)A 

contains a purely kinematic A dependence determined 
by the parton model. Letting p — ^ we obtain the result 



for b u 

^ = 2Voy''e{y) J dc^0(A^)(l + 2A^)/(^). (14b) 

Obviously, the leading-order result amounts to convo- 
luting the kinematic A dependence of the parton model 
with f{Lu). The overall factor oiy'^9{y) is not convoluted, 
since it is a trivial phase space factor, unrelated to the 
OPE. Thus, in the modified expansion we keep the factor 
y^9{y) and Eqs. are vahd (to 0(1)) over the entire 
phase space. 

Furthermore, from the first relation in Eq. © it is ap- 
parent that the twist expansion will always yield a con- 
volution of the kinematic A dependence rather than the 
rrib dependence, as originally argued in Ref. While 
a convolution of the rrib dependence is correct to leading 
order in the usual twist expansion, it introduces spurious 
subleading corrections, as has been noted before. In the 
modified expansion it already fails at leading order. 

B. Higher Twist Contributions 

The light-cone operators needed to consistently match 
all contributions of 0{K) and O(A^) are given by 



onto) = 



Here, (e • = e'^'^xnG^" (with £0123 = 1), and the 

(5-function factors are 



du!iduj2S' {to; toi, a;2)W(iX'+ + uj2)iD'^S{iD^ + uJi)j°'PLb, 

dcjidw2d(I'(5"(w; uJi,uj2, uj)b5{W+ + i^2)iD^^5{ib+ + u})W'''^ 5{ib+ + u;i)-f°'PLb, 
dtJidw2<5'(w;wi,w2)6(5(iP+ +w2)i2?*^i^'''^'5(iP+ -I-wi)7"Pl6, 
dujidu;2S'{uj] uJi,uj2)bS{iD+ + uj2){e ■ GY''5{ib+ + LUi)-f°'PLb. 

I 



(15a) 
(15b) 
(15c) 
(15d) 



S'{uj;uji,uj2) = 



6{lU — LUi) — S{lO — UJ2) 



UJi - UJ2 



(16a) 



d W2, w) = . (16b) 

OJl — UJ2 

This operator basis differs from that introduced in 
Rcf. 4] and used in previous applications 5, 7] by the dif- 
ferent Oi{llj) and the additional operator 02(1^). We note 
that this is not an artefact of our modified expansion. 
With respect to the usual twist power counting Eq. © 
both operators are formally of leading order. Neverthe- 
less, their coefficients are of at least subleading order in 
this power counting, because during the matching pro- 
cedure one effectively shifts orders from the operators to 
their coefficients by partial integration with respect to w, 
as will be illustrated later on. In turn, the B expectation 



values of 01(0;) and 02(1^) will be parametrized in terms 
of derivatives of shape functions. In the final expression 
for the spectrum these derivatives are then shifted by 
partial integration to act on the OPE coefficients. This 
ensures that in the final result the coefficients of all shape 
functions (apart from kinematic twist terms) are of usual 
twist 0(1). The same also holds for the contributions of 
©3(0;) and 'Pi{Lu) that are of usual subsubleading twist. 

The light-cone OPE of now takes the form 



+ ^ ^(2,2)) . 0^ ^ (^(1,2) ^ ^(2.2)) . 0^ 

+ (S(1'2) + ^(2,2)) . ^ 0{A^)]{UJ), (17) 
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(18a) 



where the dots denote the contraction of all Lorentz in- 
dices, and the coefficients are 

C^y\cu) = Cl^\u) = (1 - ml + 2Rl)n^ 

+ (2 - iR^ + RDA^Uo, 
C'J^l''\uj) = 2(1 - Rl){n^ + A„n„)n^ 

+ 2(1 - 3Rl + 2i?3 )77XaM, 

CZ-i^\^) = 4(1 - i?c.)^(na + A^n„)r,i^,, 
C;'(2'^)(c.) = (2(4-3i?2+2i?3K 

+ 6(2 - 2R^ + Rl)pna)n^n^ 
+ 8[(1 - R^f + 1 - 

and 

Di^Ji^) - -3(1 - R^firia + A^n„)ryi,,,, 
D'il^Huj) = -3(1 - i?^)((l + R^)n^ + 2pn^)n^n, 

- 6(1 - Rujfv±a(j,n^), 
E^^^^IH^) = -3(1 - R^fi^a + A^na)n^n,/2, 
E'Zl'^^) = 3(1 - R<^?V±al^.n,]. (18b) 

The indices enclosed in round or square brackets are com- 
pletely symmetrized or antisymmetrized, respectively. 
The ffist superscript denotes the order of the coefficient's 
term in the OPE (|17|) as it appears in the usual twist ex- 
pansion (i.e., the order of the respective shape function 
once the B expectation value is taken) , while the second 
superscript denotes the order of the coefficients's term in 
our modified expansion. 

For the above reasons, we only quote the derivatives 
of the OPE coefficients, as these are the 0{1) coeffi- 
cients which will eventually enter the energy spectrum. 
The OPE coefficients are obtained by integrating over 
uj, which increases their order in the usual twist power 

I 



counting. The constants of integration are such that each 
integral vanishes at A^j = p, that is, the kinematic 9- 
function does not contribute to the partial integrations. 

We emphasize that the OPE (fTTjl is valid to 0{A) over 
the entire phase space and to C'(A^) away from the end- 
point. When expanded into local operators, it correctly 
reproduces the full result to O(A^) J9|, as well as all local 
0{A^) terms 15j corresponding to leading and sublead- 
ing order in the usual twist expansion. 



C. Remarks on the Matching Procedure 

It is worthwhile to point out a subtlety in the matching 
procedure leading to Eq. (|17(l . The 0{A) terms contained 
in Eq. @ are 

9{A+ ~ p)C'Z''\-k+)k'^UbrPLUt, 
and can be written as a convolution in two ways, 



corresponding to the two possibilities for the matching 

C'Z^'Hu;){iD'^,6{iD++c.)}/2, CiYHc.)Or{^)- (19) 

This ambiguity has to be resolved by studying the 
one-gluon matrix element, since fc'' commutes with the 
d{k^+uj) function, while the covariant derivative and the 
d{iD+ + Lo) function do not. The gluon has momentum 
I ^ 0{A), a polarization vector e = T°-e°- , and we work 
in light-cone gauge, A+ = 0. In accordance with Eq. 
we treat Z+ exactly and expand only in L, e_, and 
ej.. To 0(A) we find 



9{A, 



-k+) - 9{A+ + ^^ 



1+) 



1+ 



e^Ubj'^PLUb + 0{A^ 



(20) 



r 



showing that we have to match onto Oi{lo). Taking the 
massless limit (i.e., p 0, R^ 0) we note that this 
is in disagreement with the results of Ref. l5t| , where the 
first possibility in Eq. (|19|1 has been chosen. Note that 
the equations of motion of heavy quark effective theory 
(HQET) cannot be used for the operator 01(11;), since the 
covariant derivative does not act directly on the heavy 
quark fields. 

A similar problem occurs in the comparison of our 
0{A^) contributions with the ones in Ref. Q. In our 
case these contributions are more complicated, requir- 
ing the two different operators C'2('jj) and 03(0;). This 



is again in contrast with Ref. 5], where only ©3(0;) ap- 
pears. In both cases the differences start at 0{A^) in 
the local expansion of the operators, and thus also in the 
final spectrum, as we will see below. 

While this paper was in the review process, studies 
of i? ^ Xu^vi based on "soft coUincar effective theory" 
(SCET) appeared which shed some light on these differ- 
ences [illilll^. The SCET-based calculations show 
that the basis introduced originally in Ref. Q is a com- 
plete basis of subleading operators, at least at tree level. 
However, in all these cases the light-cone vectors are de- 
fined based on the momentum mf,v — q, where q is the 
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total leptonic momentum. Here and in Ref. [gf a differ- 
ent choice of light-cone vectors is used, which is based on 
mtv—pi, where pe is the momentum of the charged lepton 
only. It should be possible to relate the two choices by a 
coordinate transformation, i.e., by a reparametrization. 
We shall not go into any details here, but our results 
show that for the latter choice of light-cone coordinates 
the operator basis in Ref. is incomplete. 

III. THE LEPTON ENERGY SPECTRUM 
A. Shape Functions 

In the last step we need to parametrize the B expec- 
tation values of the operators 11111 and (|15|l . To be con- 
sistent with our modified expansion we have to include 
all shape functions of leading and subleading order in 
the usual twist power counting, but also those of usual 
subsubleading twist with moments of local C(A^). 

The expansion of the QCD fields and states into HQET 
ones produces many additional operators [e.g., the Oi{uj) 
and P2{i^) of Ref s. and shape functions. However, 

these higher-order shape functions always occur in par- 
ticular combinations with those arising at leading order 
in the HQET expansion and can be suitably combined 
with them. We therefore take a different approach and 
directly parametrize the operators in QCD, which auto- 
matically combines the leading and higher-order HQET 
shape functions appropriately. 

This is in fact similar to what is used in the context 
of the local expansion, where for example the matrix ele- 
ment corresponding to the kinetic energy operator /i^ is 
also defined using the states of full QCD, and thus this 
matrix element is equal to the kinetic energy matrix el- 
ement Ai of HQET only to leading order in the 1/mf, 
expansion. 

For the leading operator we have 

{20^{iu))B - Fo{uj)v" + Ko{iu){n - (21) 



which is exact and defines the two QCD shape functions 
Fo{u!) and Kq^lo). They may be expanded into the usual 
ones of HQET, 

F„{oj) = f{oj) + ^t{oj) + 0{A^)d'{uj), (22a) 
Ko{(jj) = w/(w) + hi{oj) + 0{A^)6'{uj). (22b) 

Alternatively, we can directly perform their moment ex- 
pansions and use HQET to parametrize their moments. 
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K,{^) ^ h^hlls'i^^) + P^S"iu) + ■■■, (23b) 
where we abbreviated 

fi-fi+STa, f2 = 7-3/3 + 7-4, 
Ao = Ai -I- fi + 3(A2 + 'r2), Po^pi+^h^ 

and the pi and 71 are defined in |15| . 

We note that Fo{uj) and Ko{uj) are defined in QCD 
without any reference to the heavy quark limit. Never- 
theless, heavy quark symmetry still tells us that Ko{uj) is 
suppressed by one power of A with respect to Fo{uj). The 
normalization of Fq(u;) is exact to all orders in QCD due 
to 6-quark number conservation, while all other moments 
receive further corrections of local C'(A^) and higher. The 
leading contribution to the neglected moments is also of 
local ©(A^). 

For the higher-twist operators we find 



: -[ujFoiLu)v°' + ujKQ{uj){n - v)°']'{n - v)^' - [Fi{lo)v" + Ki{Lu){n - i;)"]'n^ - -L[{uj)r]^°''' , (24a) 
: -^[G2{io)v" + M2{Lo){n - t-)"] + \[uJ^F^{u)v'' + oj''Ko{u){n - i;)"]"(n - vY{n - vf 

+ -[luFx{lo)v°' -(- ojKx{uj){n - w)"]"2n(^(n - v)"^ -\F2(uj)v°- ^ K2{Lj){n - v)"]"n^'n'' 

+ -?7^"('^[c^ii(co)(n - vy'> + L2{u;)n'''>]" 



(:IOT{^))b = \g^{lo)v'-^^^'' - [c.2Fo(a;)]'i;"(n - vY{n -«)'' + •••, 

i:XPt^'' {uj))b = [i?4(w)(n - vY + A^4(t^)i^"]2i;[^n'^] - 77^"[^[i?4(t^)(ri - f)"' + 54(a;)n'^]]'. 



(24c) 
(24d) 



7 



These relations are again exact and define the respective 
shape functions. For the sake of completeness we give the 
full parametrization of 02(w)- In its second line we ne- 
glected all Lorentz structures whose shape functions are 
of higher order and not needed to the order we are work- 
ing. The operator Os^lo) obeys a similar parametrization 
as (!?2(^)j but with different higher-order shape functions. 

The moment expansions of the relevant shape func- 
tions are 

Li{uj) = d{uj) + —6{uj)^ , 

r< I \ 2(Ai +ri) ^pi 

- . (25) 

R^{uj) = -2(A2 + f2),5(w) - p2<5'(w) + • • ■ , 
S^[lo) = 2(A2 + t2)8[lo) + C>(A4),5'(c^) + • • • . 

The corrections to all moments shown, as well as the first 
neglected moments, are again of local ©(A"*). All other 
shape functions are of at least subsubleading twist in the 
usual power counting and do not have moments of local 



©(A^), for example 



iV4M -0(A^)5"M + -- - , 



and 



F2(w) = 0{h.^)8{uj) + 0(A4)J'H + 0(A4)5"(c^) + . . . , 
if2H = 0(A4)<5'H + -- - , 



L2M = ^<5H + o(A4)yH + ... , 

M2(tj) -C'(A'^)y'(t^) + --- . 



(26b) 



B. The Spectrum 

In order to write the spectrum in a compact way, it is 
useful to define 



r±(A) = ±(1 - 3i?^ + 2i?'') + (2 - 3i? + i?-')A, 
rf(A) = 2(l-i?3)(±i + A), 
r±(A) = ±(4-3i?2 + 2i?3) + 3(2-2i? + i?2)p^ 

r3±(A) = -3(l-i?,)[±(l + i?) + 2p]. 
The lepton energy spectrum now takes the form 



(27) 



= 2ro2/20(2/) y dt^^(A^ - p){r+(A^)^^o(^) + r;(A^)ifo(t^) + r+(A^)[c^Fo(cc>) + 2i;^i(c^)] 

+ rr(A^)[c^ifo(a') + •••]+ r+(A^)[a;2_p^(^) + •••]+ vX{^^)\lo'F^(uj) + • • •] ^^Sa) 
+ 2(1 - iRl ^ 2Rl)Li{uj) + 2(1 - R^ f{l + ^^)G2{io) 

- 3(1 - R^f[{l + A^)G3(co) - (1 - ^^)H^{io) - R^iuj) - 2Si{uj)] + ©(A^)}. 

Here, u!Ko{uj) and uj'^Fq{uj) are of usual subsubleading twist, but have a first moment of local C'(A^). The ellipsis mean 
that the same coefhcient has more shape functions, which are of the same order in the usual twist power counting, 
but have only higher-order moments, e.g., Ki{u;) or ujFi{uj). Taking the limit p — > we obtain the h ^ u result 

^ - 2Tay^e{y) J dc^0(A^){(l + 2A^)Fo(t^) + (-1 + 2A^)ifo(^) + 2(1 + A^)[LuFo{io) + 2Fi(w)] 

+ 2(-l + A^)[ujKo{uj) + ■■■]+ A[uj^Fo{uj) + ■••]- 3[oj^Fn{uj) + ■■■] (28b) 
+ 2Li{uj) + 2(1 -t- A„)G2(w) - 3[(1 + A„)G3(w) - (1 - A^)Hi{uj) - Ri{uj) - 2Si{uo)\ -f ©(A^)}. 

As for the hght-cone OPE l(T7| . Eqs. are valid to 0{K) for all lepton energies and to 0{h?) away from the 
endpoint region. They provide the correct interpolation between the two regimes of the local expansion and the usual 
twist expansion. 

Employing the moment expansions (|25(l of the various shape functions, our results reproduce the full result to local 
0{h?) Q and all local 0{K^) contributions jl5| belonging to leading and subleading order in the usual twist power 
counting. 

To compare our 5 — > u result with that of Ref. Q , we neglect all terms that are of subsubleading twist in the usual 
twist power counting and include the overall = (1 — A)^ in the power counting. Also dropping the 0{y) in front 
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and expanding -Fb(w) and Kq^lo) into HQET shape functions, we have 



^ = 2ro j dw^(A^)|(l - uj)f{Lj) + ^t{uj) - hi{cj) + 4Fi{lj) + 2G2(co) - SGsitj) + iHi{uj) 
= To j20(A) - ^[^(A) + <5'(A)] - 11A25(A) - ^ [5<5'(A) + ^^'(A)] - - y'5'(A) 



(29) 



r 



The second hne shows the expansion into local terms. 
Both expressions disagree with Ref. The differences 
arise from the new shape functions Fi{(jj) and 6*2 (w), in- 
troduced by C'i(w) and 02{^)^ and are explicit in the co- 
efficient of the /5i(5'(A) term, which is correctly contained 
in our results. In Ref. these new shape functions are 
effectively set to i^i(aj) = and G2{'-o) = G3(aj). 



IV. CONCLUSIONS 

Since ^ TibAqcD the endpoint region in i? ^ 
Xcf-Vi is affected by shape-function effects. In the present 
paper we have considered these effects to subleading or- 
der in the twist expansion. 

The usual twist expansion is valid in the endpoint re- 
gion only; however, this expansion can easily be modified 
to become valid over the full phase space, thereby yield- 
ing a smooth expression for differential rates for any value 
of the kinematical variables. We have given the relevant 
expression for the spectrum up to 0(l/m^) and to sub- 
leading order in the twist expansion. Furthermore, in a 
similar fashion as for the heavy quark expansion param- 
eters, we suggest to define shape functions using the full 
field operators and the full QCD states. 

Considering the limiting case rric — > we reveal an 
inconsistency in previous work concerning the matching 



onto subleading shape functions. It turns out that with 
our specific choice of light-cone coordinates additional 
operators are needed to obtain a complete set of sub- 
leading non-local operators. In this way, the number of 
functions needed to describe the subleading twist effects 
increases. 

The results of this paper may be useful for the es- 
timation of higher moments and for higher-order correc- 
tions to the lower moments of the lepton spectrum or the 
hadronic invariant mass spectrum in B ^ Xcii^i- Fur- 
thermore, since B Xcivg has a sensitivity to the light- 
cone distribution functions of the B meson one could also 
make an attempt to extract the shape functions from this 
process. However, for a consistent treatment one would 
need to include also radiative corrections, which could be 
considered in the framework of SCET. 
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